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Abstract 

The top-quark FCNC processes will be searched for at the CERN LHC, which are corre- 
lated with the B-meson decays. In this paper, we study the effects of top-quark anomalous 
interactions tq'y in the exclusive radiative B — >■ K*^ and B ^ pj decays. With the cur- 
rent experimental data of the branching ratios, the direct CP and the isospin asymmetries, 
bounds on the coupling kJ?^^ from B — t- K*j and kJ^^ from B ^ decays arc derived, re- 
spectively. The bound on {n^^jil from B(B — t- K*j) is generally compatible with that from 
B{B Xgj). However, the isospin asymmetry A{K*'-f) further restrict the phase of k]cRi 
and the combined bound results in the upper limit, B{t 07) < 0.21%, which is lower 
than the CDF result. For real k]^^ upper bound on B{t C7) is about of the same 
order as the 5cr discovery potential of ATLAS with an integrated luminosity of 10 fb~^. 
For B ^ decays, the NP contribution is enhanced by a large CKM factor iV^d/Vt^l, 
and the constraint on tuj coupling is rather restrictive, B{t — ^ wj) < 1.44 x 10~^. With 
refined measurements to be available at the LHCb and the future super-B factories, we 
can get close correlations between B — t- and the rare t ^ q'~f decays, which will be 
studied directly at the LHC ATLAS and CMS. 
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1 Introduction 



In the Standard Model (SM), the flavor-changing neutral current (FCNC) interactions, which 
are absent at the tree level, are highly suppressed at one loop due to the Glashow-Iliopoulos- 
Maiani (GIM) mechanism [1]. Possible new physics (NP) beyond the SM can manifest itself 
by altering the expected rates of these FCNC processes. Thus, the study of FCNC interactions 
plays an important role in testing the SM and probing NP beyond it. 

For the top quark in particular, the FCNC decays t — ?■ q'7 (where q denotes either a c- 
or a M-fiavored quark) are expected to be exceedingly rare within the SM, with branching 
ratios of order of 10"^'' [2j. Any positive signal of these decays would imply NP beyond the 
SM. These top-quark anomalous couplings could also be probed by studying the top quark 
production [3J. Present constraints on the FCNC couplings tg7 come from the following two 
experimental bounds: the branching ratio B{t — t- ^7) < 3.2% at 95% confidence level (C.L.) 
set by the CDF collaboration [4jQand the effective FCNC coupling ntu-y < 0.174 at 95% C.L. 
provided by the ZEUS collaboration These constraints will be improved significantly by 
the large top-quark sample to be available at the CERN Large Hadron Collider (LHC), which 
is expected to produce 8 x 10^ top quark pairs and 3 x 10^ single top quark annually, even 
at the initial low luminosity run (10 fb~^/year). In particular, both the ATLAS [6J and the 
CMS [7] collaborations have presented their sensitivity to these rare top-quark decays induced 
by anomalous FCNC interactions, and will be able to observe t ^ qj decays if their branching 
ratios could be enhanced to 0(10"'*) with 10 fb^^ data. 

However, when performing the study of anomalous top-quark interactions at the LHC and 

other future colliders, one should take into account constraints from some precisely measured 

low-energy processes where loops involving the top quarks play a crucial role. In this respect, 

the radiative B-meson decays induced by the FCNC transitions 6 — )■ 57 and b ^ d'y are among 

the most valuable probes of anomalous top-quark couplings [SI UHl [HI [121 [13] • For example, 

the anomalous couplings tcj mediating the rare radiative decays t ^ qj could also affect the 

-'^This upper limit corresponds to |kJ^^/A| < 1.089 TcV~^ in our convention for the effective Lagrangian 
defined by Eq. |lj). 

^This value corresponds to |k^^^/A| < 0.469 TeV~^ in our convention, and an upper limit B{t — > U7) < 0.59% 
at 95% C.L. 0. 
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radiative b ^ s'y decays through the top-quark loops [13]. On the experimental side, both the 
inclusive and the exclusive 6 — )■ 57 branching ratios are known with good accuracy (~ 5% for 
B — K*'j and ~ 7% for B — > Xs'~f), while measurements are only available for exclusive b —> dj 
channels [H]. For exclusive channels, besides the branching ratios, some other interesting ob- 
servables like the CP and isospin asymmetries have also been measured [14J. On the theoretical 
side, while the inclusive decays can be essentially calculated perturbatively, the exclusive pro- 
cesses are more complicated due to the non-perturbative strong interaction effects [15]. Besides 
some other methods [161 113 [HI [11]) the QCD factorization (QCDF) approach has provided a 
systematic framework for the treatment of exclusive radiative B- meson decays [201 [SH [23] • 
With both the experimental and the theoretical progresses, the exclusive channels are also ex- 
pected to provide important constraints on the anomalous top-quark couplings [TOl [11] and on 
various NP models 

In this paper, we shall perform a model-independent study of the effects of anomalous 
FCNC couplings tq'y in exclusive B — )■ K*'y and B ^ p'y decays, and derive the corresponding 
constraints on their strengths. Besides the branching ratios, we shall also consider the CP 
and the isospin asymmetries in these decays within the QCDF framework [201 [SB [221 [23]- 
Constraints on these anomalous couplings from the current data on these observables are then 
derived, and implications for rare t — 57 decays at the LHC are also discussed. 

Our paper is organized as follows. In Sec. |2| we introduce the effective Lagrangian describing 
the anomalous interactions tq^y, and set the convention used throughout the paper. In Sec.[3| the 
theoretical framework for exclusive B — i- V'y decays within the QCDF method, and the relevant 
formulae for rare t ^ q'y decays mediated by these anomalous interactions are presented. How 
these anomalous couplings manifest themselves in these processes is also presented here. In 
Sec. |4| we give our detailed numerical results and discussions. Our conclusions are made in 
Sec. [5j The relevant input parameters are collected in the appendix. 

2 Effective Lagrangian for anomalous tq'y couplings 

In most extensions of the SM, the new degrees of freedom that modify the ultraviolet behavior 
of the underlying theory appear only at some higher scale A. As long as we are only interested 
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in processes occurring much below this scale, we can integrate out these new degrees of freedom 
and describe the NP effects in terms of a few higher-dimensional local operators, which are built 
out of the SM fields and suppressed by inverse powers of the NP scale A [25]. A complete set of 
independent operators of dimension 5 and 6 that are consistent with the SM gauge symmetries 
could be found in Refs. [2S1 [23 EH] • 

Specific to the anomalous top-quark interactions tg7, which has a magnetic dipole structure 
as required by gauge invariance, the relevant effective Lagrangian with dimension 5 operators 
can be written in a model-independent way as [31 [281 [29] 

7 ^7 

^5 = -e 5Z ^ qR^^'tLF,, - ^ E ^ gLcr^^tiji^M^ + h.c, (1) 

q=u,c q=u,c 

where a'^'^ = i [7^,7'^]/2, ^ = (1 ± 75)5/2 are the right- and left-handed quark fields, and 
F^v = d^Ai, — duAfj_ the photon field strength tensor. Normalized to the NP scale A, the 
coefficients k^j^ and n]^^ are dimensionless couplings and in general complex. The effective 
Lagrangian given by Eq. ([T| is commonly employed in phenomenological analyses related to 
top-quark physics [21 [HI 1201 [31] . 

From the Dirac structure of the dimension 5 operators in Eq. ([T]), we can see that the tq^ ver- 
tex is induced by two independent chirality-flipped operators nig qRa^'^tLF^i, and rrit qiO'^'^tjiFf^iy, 
where the external quark mass factors must appear in order to obtain a nonzero contribution 
whenever a chirality flipping L ^ i? occurs. Due to the mass hierarchy rrit ^ rriq, the effect 
of rnqq^a^^tLFfj^y could be neglected unless the coupling kJ^^ is enhanced to be comparable to 
^ '^IqR t'y some unknown mechanism. Thus, to a good approximation, we shall focus only on 
the coupling kJ^^ in this paper [13]. 

3 Theoretical formalism 

In this section, we briefly present the theoretical framework for exclusive B — )■ V^7 decays within 
the QCDF method, and the rare t ^ qj decays mediated by anomalous tq'j interactions. For 
more details, the readers are referred to Refs. [201 [211 [221 123] and [301 [31] , respectively. 
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Table 1: Wilson coefficients at the scale = 4.45 GeV in the LL and the NLL order, using two- loop 
running for Us with the input parameters listed in the appendix. 





Ci 


C2 


Cs 


C4 






Cf 




LL 


-0.5157 


L0262 


-0.0052 


-0.0696 


0.0005 


0.0010 


-0.3179 


-0.1505 


NLL 


-0.3049 


L0082 


-0.0048 


-0.0841 


0.0003 


0.0009 


-0.3078 


-0.1692 



3.1 B — > decays within the QCDF framework 
3.1.1 The effective Hamiltonian 

In the SM, the effective Hamiltonian for radiative h D'j (with D = d,s) transitions can be 
written as 



AS^^^fi + Al-)^f,) +h.c., (2) 



where X[^^ = VqbV*j^, is the product of Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ments [32], and using the unitarity relations Xlf^ + Xi^^ + X[^^ = 0, we have 



•H^2 = C,0t + C20', + Y,C^0,, (3) 
= C,{0'i-0-,) + C2 {0^2-0^). (4) 
Here we adopt the operator basis introduced by Chetyrkin, Misiak, and Miinz (CMM) [53], 
Of = D^,{1 - 75)T>P7'^(1 - ^,)T% , Of = 57^(1 - 75)pP7ni - 75)& , 
Os = 2^7^(1 - 75)6 Yl 1^'^ ' ^5 = 2^7m7^.7M3(l - 75)& qY'Y'Y'Q , 

O, = 2^7^(1 - 1,)T% Y ^7^^^^? > ^6 = 2Dj,,j,,^,,{l - ^,)T% ^ Q^'l'^'l^'T^Q , 



^7 = - A ^^'''^(l + 75)& F^,, ^8 = - #^ Da^'^T"(l + 75)6 G;:, , (5) 



where are SU{3)c generators and rhb denotes the b-quark mass in the MS scheme. The 
corresponding Wilson coefficients at the lower scale /i = rrih can be calculated perturbatively [331 
[311 EH], and their numerical values at the leading-logarithmic (LL) and the next-to- leading- 
logarithmic (NLL) order are collected in Table [T] 
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3.2 Factorization formula for the matrix elements 

Starting from the effective Hamiltonian Eqs. ^ and (|3]), the matrix elements for B — )■ V'f 
decays can be written as [201 EI] 

(y(p',e)7(g,r/)|?^«|5(p)) = ^r«(0) |e^'^'"^r^;e>X-^ [iv*'^*) iq-p')-iv* 'P') ie*-q)] }, 

(6) 

where \V) denotes a hght vector meson state and is defined as 

(Ip^) , \K*^) , for B~ decays, 
IP />i /> y . 

-v^lpO) , , for 50 decays . 

Using this convenient parametrization, all the dynamical information is encoded in the form 
factor 7^*^(0), which constitutes a big challenge to be determined precisely [201 Ell 122] . 

In the QCDF formalism, the form factor 7^*^ (0) could be computed in terms of heavy-to- 
light transition form factors and hadron light-cone distribution amplitudes (LCDAs) at the 
leading power in l/m^ expansion. Explicitly, it takes the following factorization formula [20] 

r«(0)=Ti(0)C« + ^^$^ / fduMu)Tf^{u,uj), (8) 

iVc niB ^ J ^ Jo 

where J'b and $ b,± denote the B-meson decay constant and LCDAs, and fy and (j)± the corre- 
sponding quantities of the light vector meson, respectively. The first term, expressed in terms of 
the tensor form factor Ti{0), corresponds to the vertex corrections where the spectator quark 
in the B-meson does not participate in the hard process, whereas the second term incorpo- 
rates the hard-scattering contributions where the spectator quark is involved. Thus, the first 
and the second term are usually called the "form factor" and the "spectator scattering" term, 
respectively. 

The hard-scattering kernels C^'' and in Eq. (8) can be calculated perturbatively and, 
up to next-to-leading order (NLO), have the following expansions, respectively [20] 

C« = Cf'''> + ^ d^''^ + ..., (9) 

4:71 

Tfl^ {u, CO) = Tf;2 (u, ^) + ^ T[]i (n, c.) + . . . , (10) 
where the QCD coupling should be evaluated at the scale /if, ~ rrih in Eq. ^ and at 



— (mfeAqcD)^^^ in Eq- (10), corresponding to the typical virtualities in the two terms, 
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respectively. Explicit expressions for the coefficients C^^'^\ C^'^^ and could be 

found in Ref. EOl. 



For B —7- V'y decays, besides the leading-power contributions given by Eqs. (|9]) and (10), 
it is also well-known that some very specific power corrections, such as the weak annihilation 
and exchange amplitudes in B ^ that are enhanced by the large Wilson coefficient C2 ~ 1, 
are often numerically important [20l |2T| [22]. Furthermore, the annihilation topologies have 
been shown to provide the main source for the isospin asymmetry in B K*'y [23 j. Although 
being power-suppressed by l/rrib, these contributions are still computable within the QCDF 
framework. Thus, in this paper we have also included these isospin-breaking power corrections, 
which are denoted by A7^*''|ann and A7^*''|hsa for the weak annihilation and the hard spectator 
scattering, respectively. Their explicit expressions could be found in Refs. [20| [36]. 

For the isospin-breaking power corrections, it should be noted that an endpoint divergence 
is encountered in the matrix element of chromo-magnetic dipole operator Oq, belonging to the 
term AT^^^sa. Following the treatment adopted in Refs. [221 [3S], we regulate this singularity 
with an ad hoc cutoff 

du^ {l + pe''^) / du, (11) 



'0 ^0 

and take A/j ~ 0.5 GeV, together with < p < 1 and < < 27r, to give a conservative 
estimation of the theoretical uncertainty related to this power correction. 

3.2.1 Observables in i? — > V'y decays 

When discussing observables in i? — )■ decays, it is more convenient to express the decay 
amplitudes in terms of a new quantity defined by [20] 

^7 - ^7 +••• , (12) 

where Cf = - C^/?, - ACJ^ - 2OC5/3 - SOCg/Q, defined in the MS scheme with fully 
anti-commuting 75, is the effective Wilson coefficient [33], and the ellipses denote the 0{as) 
and the sub-leading power corrections discussed in the last subsection. 

In terms of the quantity Cj \ the decay rate for B — t- V'y decays can then be written as [20] 

r(5 ^ V^7) = ^ -I 5 f 1 - ^ rnl UOf W^C? + M^^Ci^^l' , (13) 



with S = 1/2 for p°, and 5 = 1 for the other hght vector mesons. In the SM, the CP- 
conjugated mode T{B — t- V7) follows from Eq. (13) with the replacement A^-^^ — )■ A^-'^''*. For 



h ^ s transitions, the dominant term is Aj^'^C^*'*, since the contributions proportional to are 
doubly Cabibbo-suppressed; whereas for h ^ d transitions, where X^'^ is of the same order as 
Xf'\ the interference term is non-negligible and can be the source of interesting CP- violating 
and isospin-breaking effects. 



With the decay rate Eq. (13) at hand, the interesting observables in i? — V^7 decays can 



be defined as follows [201 |2ll 122] 

• the CP-averaged branching ratio 



with tb being the i?-meson lifetime, 
the direct CP asymmetry 



A n/ ^ r(^ -> v^) - T{B ^ v^) 



the isospin asymmetry 



* " r(B»-»ir*»7) + r(B+-»A'*+7)' ^ ' 

where all decay rates are assumed to be CP-averaged. 

There are five observables for B — )■ K*'y and B ^ decays, respectively, i.e., two CP- 
averaged branching ratios, two direct CP asymmetries, and one isospin asymmetry. Once being 
measured precisely, they could be used to test the SM and to probe various NP beyond it [101 
[TTl 121]. Especially, the two isospin asymmetries are expected to provide useful information 
complementary to the corresponding inclusive decay modes 



3.2.2 Anomalous igj coupling effects on i? — V'j decays 

The anomalous tq'~f interactions affect b — t- -D7 transitions through the two Feynman diagrams 



depicted in Figs. 1(a) and 1(b) It is interesting to note that the associated CKM factors in 



(a) 



(b) 



(c) 



Figure 1: Feynman diagrams for b — )• s{d)'y transitions in the unitary gauge, (a) and (b) are mediated 
by the anomalous tq^ interactions defined by Eq. ([T]), while (c) represents a sample LO penguin diagram 
in the SM. 



Figs. 



la 



and 



1(b) are VtbV*j^ and K^^V^^, respectively. Since |VjfeVg*£,| ^ iKb^ti)! ^"^^ both 



q = u and c quarks, the contribution of Fig. 1(a) would be much larger than that of Fig. 1(b) 



Furthermore, given the strengths of the couplings twy and tcj comparable, the contribution 



of Fig. 1(a) is still dominated only by one coupling, tcy for 6 — 57 and twy for 6 — ^7 



respectively, because of the relations \Vcs\ > \Vus\ and \Vud\ > \Vcd\- Hence we shall only 



consider the contribution of Fig. 1(a) with only one anomalous coupling. 



From the Feynman diagram Fig. 1(a), it is also observed that the large CKM factors 
l^bKsl ~ 1 and |V(feVud| ~ 1, compared to the corresponding SM case jVibVisI ~ (9(A^) and 
|14bVjd| ~ (9(A^), make the transitions 6 — 57 and 6 — )■ ^7 to be very sensitive to the strengths 
of anomalous couplings tcy and twy, respectively. Constraint on the coupling tcy from the 
precisely measured inclusive decay B — t- Xg'y has been studied in detail by two of us 



The calculation of Fig. |l(a) can be most conveniently carried out in the unitary gauge where 
the pseudo-Goldstone components of the SM Higgs doublet are absent, which has been calcu- 
lated in 't Hooft-Feynman gauge in Ref. pS]- It is noted that in the unitary gauge, the g"^^ part 
of the ly-boson propagator gives an ultraviolet-finite contribution, while the contribution from 
the q'^q^ /rn^ part is ultraviolet-divergent. Following the treatment adopted by Grzadkowski 
and Misiak [12], we shall apply the modified minimal subtraction (MS) scheme to absorb the 
divergences into some counterterms served by other dimension-six operators [2S1I2Z1, the MS- 
renormalized Wilson coefficients of which are assumed to be negligible in comparison with the 
ones considered here. With such a prescription, logarithms In are present in the matching 
coefficients. 
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Normalized to the effective Hamiltonian Eq. ^ and the operator basis Eq. ([5]), the contri- 
butions of anomalous tq'y interactions to 5 — )■ V'~f decays would result in the deviation 



■yNP/ 



(18) 



where {xg = mq{^,wY /m^) 



ts 



In 



4 ^ 2(Xr 



+ 



2{xc - iy{xc - Xt) 



InXc + 



ri 



2{xt - lY{xt - Xc 



In Xt 



(19) 



for b ^ S'-f transition, and 



w 



K 



tuR 



ud 



td 



ln^-- + 



1 



+ 



mw 4 ' 2{xu - l){xt - 1) 

3 

lnx„ + 



X 



xt 



In Xt 



2(x„ - l)2(a;„ - Xt) ' 2(xt - l)2(xt - x„) J 

for b ^ dj transition, respectively. It is noted that the NP contribution C^^i^fiw) is suppressed 
by a mass factor rrit/A, but enhanced by a CKM factor V*j^/V^*^. Since the NP contribution 
does not bring about any new operators, the renormalization group evolution of the Wilson 
coefficients from the scale fiw down to fib is just the same as that in the SM. 

As a final remark, we also find that the operator qRcr^'^tLFf^^ in the effective Lagrangian 
Eq. ([l]) contributes to 6 — > D'y transitions only through a term ra£,Da^y{l — 75)^. Neglecting 
the light quark mass m^), as done in the SM, the effect of qRCT^'^tLF^y could be therefore safely 
neglected, which supports the remarks on this operator made in Sec. [2] 



3.3 Rare t ^ decays mediated by anomalous tq^ interactions 

Since t — )■ bW is the dominant top-quark decay mode, the branching ratios of radiative t ^ 
decays are usually defined as 



B{t ^ qj) 



T{t ^ g7) 



(21) 



T{t bW) ■ 

As the SM predictions for r(t — )■ 57) are exceedingly small [2], we need only consider t — )■ 57 
decays mediated by the anomalous tg7 interactions, which have been recently calculated at the 
NLO in Refs. [30l [31]. The final result for the decay widths T(t — > qj) reads [301 EI] 



r(t^g7) = ro(t^g7) 1 



2a^ 

97T 



-31n(^' 

ml 



27r^ + 8 



(22) 
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where ro(t — )■ 57) = aeml \k1^^/K\'^, is the leading order (LO) decay width. 

The decay width of the dominant top-quark decay mode t — )■ hW at the LO and the NLO 
could be found in Ref. and is given below 



where ro(t ^ hW) = ^^\Vtb\'^l3^{3-2l3l.), is the LO decay width and (3w = {l-ml^/mfy/^, 



is the velocity of the ly-boson in the top-quark rest frame. 

4 Numerical results and discussions 

4.1 Numerical results for the observables in 5 — > decays 

With the theoretical framework discussed above and the input parameters collected in the 
appendix, we first present the SM predictions for the CP-averaged branching ratios, the direct 
CP asymmetries, and the isospin asymmetry in B K*'~f and B ^ p'j decays, respectively, 
which are listed in Table |2} The theoretical uncertainties are obtained by varying the input 
parameters within their respective ranges. The experimental data are taken from the Heavy 
Flavor Average Group pH] . 

From Table [2| we can see that, with their respective uncertainties taken into account, the SM 
predictions for the branching ratios are in good agreement with the experimental measurements. 
The predicted signs of the direct CP asymmetries in B ^ K*'~f decays are opposite to the 
experimental measurements, but with large uncertainties. It is noted that, within the QCDF 
formalism, the main theoretical uncertainty for the branching ratio comes from the tensor form 
factor, whereas for the direct CP asymmetry, the largest theoretical uncertainty is due to the 
residual renormalization-scale dependence (for B — )■ p7) and the first Gegenbauer moment 
tti (for B — 7- K*'~f). Since the isospin asymmetry in these decays is a power-suppressed effect, 
its calculation is less certain; however, this observable can usually provide useful information 
about NP parameter spaces complementary to the inclusive mode B — )■ [21] 

As the current experimental data and the theoretical predictions for these observables still 



r(t ^ bW) = To{t ->bW){l + 



3-2/3i 



9-4/3; 




(23) 
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Table 2: Observables in S — )• K*^ and B ^ decays. The branching ratios are given in unit of 
10~®, while the CP and isospin asymmetries are given in unit of 10~^. 



Observables 


Exp. data ^ 


SM prediction 


B{B+ K*+-f) 


42.1 ± 1.8 


44.60;j?J« 




43.3 ± 1.5 


45.99tJ?i° 




18 ±29 


-n 1 5+°-2^ 

o.iO_Q 22 


Acp{K*'l) 


-16 ±23 


U.OU_Q 29 




5.2 ±2.6 


c, 1 +2.83 


B{B+ p+7) 


0.98l°J^ 






0.86l°:J^ 


o.8ol°:i 


•Acp{p^l) 


-11 ±33 


^•"^^-3.85 


^cp(p°7) 




-8.9412-26 


A(p7) 




-8.401^:^1 



have large uncertainties, in the following numerical analyses, we shall consider them with la 
theoretical and 2a experimental uncertainties, respectively. 

4.2 The anomalous coupling k:^^^ in exclusive B — > A'*7 decays 

For B — )■ i^*7 decays, the main contribution is due to the anomalous coupling kJ^^. With the 
notation kJ^^ = e*^t<:«, we get numerically 

C;f_,^(/i,) = -0.3179 ± 2.3985 e*(-i78.96°+e- ^) ^ (24) 

at the LL approximation, where the NP scale A is given in unit of TeV. This indicates that, for 
a given value |kJ^^|/A, the NP contribution is constructive to the SM one in the region 9]^j^ ~ 0°, 
whereas in the regions ^^cr — ±180°, the interference between them becomes destructive. 

The upper bounds on the anomalous coupling |/tJ^j:j/A| as a function of 9]^^, constrained 
by the five observables in i? — )■ K*^ decays, are shown in Fig. |2} The light and the dark blue 
region are obtained with la theoretical and 2a experimental uncertainties, and with only 2a 
experimental uncertainty, respectively. The dashed curves are obtained with central values. 
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-150 -100 -50 50 100 150 



(a) B{B+ K*+j) 



(b) B{B" K*^-f) 



0.4 I 

> 0,3 f 
H 

^ 0,21- 



01 

0,0 F 



-150 -100 -50 50 

Cr(°) 



(c) Acp{K*+j) 



100 150 



> 
H 

0.2 - 

0,1 

0,0 



-150 -100 -50 50 



(d) Acp{K*"^) 




100 150 




-50 50 

(e) A(i^*7) 





-50 50 



(f ) Combined constraints 



Figure 2: The upper bounds on the anomalous coupUng |kJ^^/A| as a function of Oj^j^, constrained 
by the five observables in i? — )• K*j decays. The Hght and the dark blue region are obtained with la 
theoretical and 2a experimental uncertainties, and with only 2a experimental uncertainty, respectively. 
The dashed curves are obtained with central values of the input parameters. 
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From Figs. 2(a) and 2(b) , we can see that the constraints from the two branching ratios are 
quite similar to, but shghtly looser than the one from B{B — )■ Xg'y) [12]. This is due to the 
fact that both the exclusive and the inclusive branching ratio are, at the LO approximation, 
proportional to the same Wilson coefficient jCyp. The large theoretical uncertainty is also clear 
by comparing the light blue region with the dashed curve. In addition, the region 6^^^ ~ 0, 
where the NP contribution is constructive to the SM one, gives the most stringent upper 
bound on the magnitude |kJ^^/A|. There are, however, two allowed solutions in the regions 
— ±180°, the larger one corresponding to the case in which the sign of C^^ is flipped. 

Due to the large theoretical and experimental uncertainties, the two direct CP asymmetries 



currently provide no valuable information about the anomalous coupling, as shown in Figs. 2(c) 



and 2(d) However, once measured and computed precisely, as shown by the dashed curves, 
these observables are also expected to provide very useful constraints. 

Although being a power-suppressed effect, the isospin asymmetry, being proportional to Cf ^ 
to first order, can provide complementary constraints on the anomalous top-quark coupling, as 



shown in Fig. 2(e) Especially, the large sign-fiipped solution allowed by the two branching ratios 
can be, to some extent, reduced. Since we have included 2a uncertainty of the experimental 
data, the constraint from this observable is quite loose, particularly in the region 6^^^ ^ 0. 
Combining all the five constraints from B — K*^ decays, we present our final results in 



Fig. 2(f) We can see that there is actually no constraint on the phase 6tcR, while the upper 
bound on the strength |kJ^^/A| is more stringent than obtained from the inclusive B — )■ Xs'j 
decay. The implication for radiative t — )■ 07 decay will be discussed later. 



4.3 The anomalous coupling kJ^j^ in exclusive B ^ decays 

For B ^ p'j decays, the main contribution is due to the anomalous coupling k^^^, and we get 
numerically 

C',f_,,^{fi,) = -0.3179 + 11.3477 e^(-2^-^8°+e7„«) ^ (25) 



at the LL approximation. In contrast to Eq. (24), the destructive interference between the 



NP and the SM contribution occurs in the regions ^J^^ ^ 22°. The NP contribution is also 
accompanied with a larger coefficient, which is due to the larger CKM factor |\4d/V(d|. 

As done in the case of i? — ?■ K*^ decays, we show in Fig. [3] the upper bounds on the 
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(e) Combined constraints 

Figure 3: The upper bounds on anomalous coupling {kJ^j^/AI as a function of ^7tj_R' constrained from 
the four observables in i? — )• p7 decays. The other captions are the same as in Fig. [2} 

anomalous coupling |/t7„^/A| as a function of 0]^j^, constrained by the four observables in 



B ^ decays. Due to the numerical difference between Eqs. (24) and (25), the behavior of 
the constraints on the NP parameters is quite different. 
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As shown in Fig. 3(b) , constraints from the branching ratio — i- p°7) exhibit the 
destructive (near the region 6].^^ ~ 22°) and constructive (near the region 6].^^ ^ —158°) 
features between the NP and the SM contribution. The loose bound in the destructive region 
corresponds to the case in which the sign of Cj^ is flipped. On the contrary, since the SM 
prediction is much larger than the experimental data (see Table [2]), the constructive region is 



already excluded by the branching ratio B{B^ — )■ p^7), as shown in Fig. 3(a) Thus, the phase 
of the anomalous coupling 6]^p^ is constrained into a small region by this observable. 

Similar to the case in i? — )• K*'y decays, the direct CP asymmetries, with the current 
experimental measurements, could not give any constraints on the anomalous coupling kJ^^, 



which is shown in Fig.|3(c) With respect to the isospin asymmetry A(p7), as shown in Fig. 3(d) 



it is found that the upper bounds on |fi;J^^/A| are more restrictive in the region 6'/^^ G [0, 180°], 
while quite loose in the region 9]^j^ G [—180°, 0]. It is, however, expected to provide interesting 
constraints once the data becomes more precise. 



The final combined constraints are shown in Fig. |3(e) Although being excluded in the 
region 9]^j^ G [0, 180°], the large sign-flipped solution is still left in the region 9]^j^ G [—180°, 0]. 
As already discussed, being enhanced by the large CKM factor V*^/V^*^, the anomalous coupling 
'^tuR is expected to be severely constrained by these decay modes. However, due to the large 
theoretical and experimental uncertainties, the obtained upper bound on is at the same 

order as obtained from B K*'y decays. 



4.4 Combined upper bounds on tqj and implications for B{t — > q^) 

Taking 6]^^ = 0° and ±180° as benchmarks and setting A = 1 TeV, we summarize our numerical 
constraints on the strength \n]qji\ in Tables [s] and |4| where SI and S2 correspond to the two 
cases in which the sign of C^^ is not flipped and flipped, respectively. For comparisons, the 
constraint from B — )• Xgj [IB|, and the upper bounds at 95% C.L., [kJcrI < 1-089 from CDF [1] 
and < 0.469 from ZEUS [5], are also given. 

From Tables [3] and |4| we can see that, the strengths |kJc_rI I'^IurI constrained by these 
exclusive decay modes are both lower than the current experimental limits. As shown in the 
last column of Table [3| the large sign-flipped solution allowed by the two branching ratios does 
not survive the constraint from the isospin asymmetry A{K*'j). It is also noted that, for a real 
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Table 3: Constraints on from B — )• K*^ decays, with fixed and A = 1 TeV. SI and 

S2 correspond to the two cases in which the sign of Cy^ is not flipped and flipped, respectively. For 
comparisons, the constraint from B — )• Xg'y [l3] and the upper limit from CDF [4J are also listed. The 
corresponding upper bounds for B{t — )• cy) for each case are listed in the last row. 





Oln = 0° 


= ±180° SI 


elj, = ±180° S2 




[0, 0.018] 


[0, 0.032] 


[0.311, 0.361] 




[0, 0.017] 


[0, 0.032] 


[0.329, 0.379] 






[0, 0.096] 


[0.175, 0.177] 


Combined 


[0, 0.017] 


[0, 0.032] 







[0, 0.016] 


[0, 0.019] 


[0.45, 0.48] 


CDF bounds [U 


[0, 1.089] 


[0, 1.089] 


[0, 1.089] 


B{t^ C7) 


< 7.79 X 10-6 


< 2.76 X 10-5 






Table 4: Constraints on from B ^ p'j decays, with fixed 0^^^ and A = 1 TeV. For a 

comparison, the upper limit from ZEUS ^ is also shown. The other captions are the same as in 
Table [S 





^tuR 


= 0° SI 


OIr = 0° S2 


= ±180° 


B{B+ p+7) 


[0, 


0.020] 


[0.051, 0.071] 







[0, 


0.011] 


[0.060, 0.081] 


[0, 0.010] 


Combined 


[0, 


0.011] 


[0.060, 0.071] 





ZEUS bounds ^ 


[0, 


0.469] 


[0, 0.469] 


[0, 0.469] 


B{t wj) 


< 3.26 X 10-6 


[9.71 X 10-^ 1.36 X 10"^] 






coupling kJc/j, the inclusive B — t- Xs7 decay provides more restrictive bounds than the exclusive 
B K*'y decays. 

Finally, we collect in Table |5] the combined upper bounds for the anomalous couplings tqj, 
and the corresponding predictions for B(t — )■ ^7). The dependence of B(t -> ^7) on the phase 
9]^j^ is shown in Fig. |4| where for comparisons, the 95% C.L. upper limit B{t q^y) < 3.2% 
from CDF [3] and B{t M7) < 0.59% from ZEUS [5J, as well as the ATLAS sensitivity with a 
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Table 5: Combined upper bounds for the anomalous couplings tq^ from exclusive radiative B-meson 
decays, and the corresponding predictions for B{t — >• 57). For comparisons, we also list the upper 
bound from B Xsj pj], CDF |31 and ZEUS [5], respectively. 





B K*-f 


5 ^ X,7 [T3] 


5 — > P7 


^tcR [°] 






din n 


[-71.6,100.8] 




< 0.277 


< 0.48 




< 0.073 


B{t C7) 


< 0.21% 


< 0.63% 




< 1.44 X 10-^ 


B{t C7) CDF j4j 


< 3.2% 




B{t^ wy) ZEUS |5] 


< 0.59% 
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?; 0.001 - 

T 

§ 10-" t- 
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(a) B{t C7) 



(b) B{t uj) 



Figure 4: The upper bound on B{t qj) as a function of 6*7^^. Both the CDF [4| and the ZEUS [5] 
upper limit are given at 95% C.L.. The ATLAS sensitivity [6], with a 5a significance at an integrated 
luminosity of 10 fb~^, is also shown. 

5a significance [H], B(t — )■ ^7) < 9.4 x 10~^, are also given. 

From Table [sj we can see that the combined upper bound on obtained from B — t- K*'y 
decays is a bit smaller than that from the inclusive B — )■ Xg'j decay. As mentioned already, this 
is mainly due to the constraint from the isospin asymmetry, which makes the large sign-flipped 
solution reduced. As shown in Fig. |4| the predicted upper limits on B{t — )■ 07) and B{t — wy) 
are lower than the ones from the CDF [4J and the ZEUS [5J collaboration, respectively. They 
are, however, of the same order as the 5a discovery potential of ATLAS [6j, with an integrated 
luminosity of 10 fb~^. 



18 



5 Conclusions 



In this paper, within the QCD factorization formahsm, we have studied the effects of anomalous 
top-quark FCNC interactions tq^ in exclusive radiative B K*j and B ^ pj decays. Among 
the two dimension-5 operators, QRcr^'^tLFi^,^ and ^iCr'^^iflF^j,, only the second one is found to give 
a main contribution to these decays. With the current experimental data of the branching ratios, 
the direct CP and the isospin asymmetries, bounds on the couplings i^4qW which determines 
the strength of this operator, are derived. Our main conclusions are summarized as follows. 

For B K*j decays, the main contribution is due to the coupling k^^^. The combined 
constraints from the two branching ratios and the isospin asymmetry exclude the region where 
the sign of Cf^ is flipped due to the NP contribution. As a result, the upper bound on {kJ^^,/ A\ 
is stronger than that obtained from the inclusive B — > Xgj decay. The corresponding upper 
limit, B{t — )> C7) < 0.21%, is lower than the CDF result. However, for real k.]cr, the bound 
is rather restrictive, and the corresponding upper limit on B{t ^ cj) is generally of the same 
order as the 5a discovery potential of ATLAS with an integrated luminosity of 10 fb~^. 

For B ^ fTf decays, the main contribution is, on the other hand, due to the coupling nj^^ 
and enhanced by a large CKM factor \Vud/ytd\- Due to the large theoretical and experimental 
uncertainties, constraints on the anomalous coupling are not as strong as expected. However, 
with the current data and our SM prediction of B{B^ — >■ ^"'"7), between which there is a large 
difference, the phase of the anomalous coupling O^^j^ gets, to some extent, constrained. The 
obtained upper limit, B{t — > ^7) < 1.44 x 10~^, is stronger than the ZEUS result. For most of 
the constrained parameter space, the upper limit on B{t — )■ ^7) is about the same order as the 
5(7 discovery potential of ATLAS with an integrated luminosity of 10 fb~^. 

As a result, there is an interesting interplay between exclusive radiative B-meson and rare 
i — > Q'7 decays. With refined measurements from the LHCb and the future super-B factories, 
we can get close correlation between the 6 — > S7 decays and the rare t ^ q'y decays, which will 
be studied directly at the LHC CMS and ATLAS. 
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Appendix: Theoretical input parameters 

In this appendix, we collect all the relevant input parameters, when calculating the observables 
in i? — 7- V'-f decays and the branching ratios of t — g7 decays mediated by the anomalous tq^ 
interactions. 

The basic SM parameters 

First, we need some basic SM parameters, which are, if not stated otherwise, taken from the 
Particle Data Group [8] 

asimz) = 0.1187 ± 0.0007, a = 1/137.036, Gp = 1.16637 x 10"^ GeV"^ 
sin^Ow = 0.23146, mw = 80.399 GeV, mz = 91.1876 GeV, rrit = 173.3 ± 1.1 GeV [3B], 
mB+ = 5279.17 MeV, uibo = 5279.50 MeV, mK^+ = 891.66 MeV, m^^o = 895.94 MeV, 
mp+ = nipo = 775.49 MeV, tb+ = 1.638 ps, tbo = 1.525 ps, (26) 

where rrit is the top-quark pole mass. We use two-loop running for as throughout this paper. 
The CKM matrix elements 

For the CKM matrix elements, we adopt the Wolfenstein parametrization and choose the 
four parameters A, A, p and rj as fitted by the CKMfitter group [30] 

A = 0.812;n^?, A = 0.22543l°:S , P = 0.144l°:°i , 7^ = 0.342;°;°}^ (27) 
with p = p{l — ^) and f] = i] {1 — 
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Quark masses and meson parameters 

For the running quark masses in the MS scheme, we take [H] 

fribimb) = 4.164 ± 0.023 GeV, mc(mc) = 1.273 ± 0.006 GeV. (28) 

To get the corresponding pole and running quark masses at different scales, we use the NLO 
MS-on-shell conversion and running formulae collected, for example, in Ref. |12] . 



It should be noted that the quantity C^^\ defined by Eq. (12), depends on the 6-quark mass 
renormalization scheme. Following Ref. [20], in this paper we use the potential-subtracted (PS) 
scheme [13]. Converting mb{fnb) — 4.164 GeV to the PS mass mps(2 GeV), we get numerically 
mpg(2 GeV) ~ 4.45 GeV. To give an estimate of renormalization scale uncertainty, the two 
renormalization scales Hb and Hhc are chosen, respectively, as 

l^b = 4.45^^;^^ GeV, i^hc = 1-5 ± 0.6 GeV. (29) 

When discussing exclusive B — )■ V'j decays, we need some parameters related to the involved 
mesons, such as the transition form factors, the decay constants, as well as the Gegenbauer 
moments of meson LCDAs. Most of these parameters are not directly known from experiment, 
and have to be determined by some nonperturbative methods like the QCD sum rule and lattice 
QCD. A summary of these parameters is listed below 

Tf^^*(0) = 0.31 ±0.04 [18], Tf^''(0) = 0.27 ±0.04 [18], /^^ = 192.8 ± 9.9 MeV [Sj, 

Aij,+ (1.5 GeV) = 485± 115 MeV [20j, //^* = 220 ± 5 MeV [18], /p = 216 ± 3 MeV [18], 

/^,(1 GeV) = 185 ± 10 MeV [IS], fp{l GeV) = 165 ± 9 MeV ^3], 

ai(K*)x,||(l GeV) = -0.04 ±0.03 [18], a2(K*)±,||(l GeV) = 0.15 ±0.10 [H], 

a2(p)±,||(l GeV) = 0.15 ± 0.07 [E], (30) 

where the tensor form factors Ti evaluated at = 0, the decay constants and Gegenbauer 
moments of the light mesons are obtained from QCD light-cone sum rules (LCSRs) [TH]- The 
quantity Xb,+ is defined as the first inverse moment of the B-meson LCDA [:20j, and is the 
B- meson decay constant. 
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